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Abstract 



A new integrable class of Davey-Stewartson type systems of nonlinear partial diffe- 
rential equations (NPDEs) in 2 + 1 dimensions is derived from the matrix Kadomtsev- 
Petviashvili equation by means of an asymptotically exact nonlinear reduction method 
based on Fourier expansion and spatio-temporal rescaling. The integrability by the 
inverse scattering method is explicitly demonstrated, by applying the reduction tech- 
nique also to the Lax pair of the starting matrix equation and thereby obtaining the 
Lax pair for the new class of systems of equations. The characteristics of the reduc- 
tion method suggest that the new systems are likely to be of applicative relevance. 
A reduction to a system of two interacting complex fields is briefly described. 



1 Introduction 

> 

New classes of evolution nonlinear partial differential equations (NPDEs) integrable by 
the inverse scattering method (S-integrable) have been found in the last years. These 
equations are known to be applicable to various branches of physics such as fluid dynamics, 
nonlinear optics, condensed matter physics and so on. The most famous examples are the 
Korteweg-de Vries and the nonlinear Schrodinger equations in 1 + 1 dimensions and the 
Kadomtsev-Petviashvili and the Davey-Stewartson equations in 2 + 1 dimensions [1]. 

A simple explanation of this coincidence (integrability and applicative relevance) is 
based on the observation that very large classes of evolution NPDEs in 1 + 1 and 2 + 1 
j> \ dimensions, with a dispersive linear part, can be reduced, by a limiting procedure involving 

the wave modulation induced by weak nonlinear effects, to a very limited number of 
"universal" evolution NPDEs. Moreover, the same model equations obtained in this way 
appear in many applicative situations (for instance in plasma physics, nonlinear optics, 
hydrodynamics, etc.), where weakly nonlinear effects are important [2-5]. 

The reduction method preserves integrability and therefore the model equations are 
likely to be integrable. For example, it is sufficient that the very large class of equations 
from which they are obtainable contains just one S-integrable equation, provided the lim- 
iting procedure preserves integrability, so that the property of S-integrability is inherited 
through this limiting technique. Obviously, the last statement about the integrability is 
based on heuristic considerations and could not be characterized as a rigorous theorem. 
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No precise definition of integrability is available for evolution NPDEs, there being much 
difference between finding the general solution of a NPDE or solving an initial- value prob- 
lem with given input data and boundary conditions. It would be possible to derive the 
spectral transform of the Davey-Stewartson equation from the spectral transform of the 
Kadomtsev-Petviashvili equation. 

Thus this approach, besides explaining why certain model equations are integrable and 
applicable, provides a powerful tool to investigate the relation among different integrable 
equations, to test the integrability of nonlinear evolution PDEs and, most importantly, to 
identify integrable evolution equations that are likely to be of applicative relevance. 

In previous papers, we applied this method to certain integrable equations in 2 + 1 
dimensions. The most interesting results are that the Davey-Stewarston equation [6-7] 
is the typical model equation in 2 + 1 dimensions, while new integrable NPDEs can be 
obtained together with their Lax pair [8—11]. Moreover, we used the reduction method to 
derive two equations of applicative relevance in plasma physics [12-13]. 

The basic idea of the reduction method is to consider a nonlinear evolution PDE whose 
linear part is dispersive; as it is well known the linear evolution is most appropriately 
described in terms of Fourier modes and each Fourier mode evolves with constant am- 
plitude and an associated group velocity, that represents the speed with which a wave 
packet peaked at that Fourier mode would move in configuration space. To evaluate the 
weak nonlinear effects it is convenient to consider a specific Fourier mode and follow it by 
going over to a frame of reference that moves with its group velocity. The weak nonlinear 
effects give rise to a modulation of the amplitude of that Fourier mode (that would remain 
constant in the absence of nonlinear effects). The modulation is best described in terms 
of rescaled "coarse-grained" and "slow" variables, that display the weak nonlinear effects 
on larger space and time scales; indeed, the first step of the reduction method is to use a 
moving frame of reference with the introduction of the slow variables: 

£ = £ P(x - V ± t), r, = e p {y - V 2 t), r = e% 

p>0, q>0, (1.1) 

where V\ = V±(Ki, K2), V2 = V2 (if 1,^2) are the components of the group velocity 
V_ (K_) = (V\(Kx, K2), V2(Ki, K2)) of the linearized equation, i.e. of the equation obtained 
by neglecting all the nonlinear terms, and e is a "small" expansion parameter. 

It is thereby seen that the function that represents the amplitude modulation satisfies, 
in terms of the rescaled, slow, variables, evolution equations having a universal character; 
since the coarse-grained nature of the new variables implies that only certain general 
features of the nonlinear interaction are important. 

In this paper we expose an interesting extension of this approach and consider the 
matrix Kadomtsev-Petviashvili equation [14-15] 

U t + U xxx -W y + iV3 [W, U] - 3{U, U x } = 0, 

W x = U y , (1.2) 

where [A,B] = AB — BA, U = U(x,y,t), W = W(x,y,t) are N <g) N complex matrices 
and the subscripts denote partial differentiation. 
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By applying the reduction method, a new class of integrable matrix systems of evolution 
NPDEs depending on a real parameter A is obtained 

i^ T + L^- A[tf,A] - + V3 {A, ^} - {$,^ 2 } = 0, 

- L$ - A [$, A] - [17, $] - ^3 {A, $} + {*, $ 2 } = 0, 
(3 - A 2 ) A 5 + 2AA,, - n v - V3{*, + [V, + A [$, *] € = 0, 

A„ = % (1.3) 

where {^4, 73} = >1_B + Bvl, the linear differential operator L is given by 

i = -(3+^)^ + 2A^-^, (1.4) 

and * = ^(£,77, r), $ = $(£,77,7-), A = A(£,?7,t) and O = fi(£, 77, r)are N ® N complex 
matrices. 

The paper is organized as follows. In the next section we apply the reduction method 
to the starting equation ( |1.2| ) and obtain the new system of matrix equations (|1.3|)-(p~i|). 
Moreover, we reduce the matrix system of equations to a new integrable two-component 
complex fields system of nonlinear equations, which, in the one-component case, reduces 
to the standard Davey-Stewartson equation. In Section 3 we discuss in some detail how 
the reduction method can be applied to the Lax pair of the equation (|1.2|) and we derive 
the Lax pair of the system of matrix equations ( |1.3| )-( |L4| ). Finally in the last section we 
recapitulate the most important results and indicate some possible extensions. 

2 A new integrable matrix system in 2 + 1 dimensions 



Lemma 2.1. The linear dispersive part of the starting equation admits as a solution 
a Fourier mode, with a group velocity V_(K_) = (Vi(Ki, K%), V%(Kx, -K^))? 



V X {K X ,K 2 ) = -ml + |f, V 2 {K U K 2 ) = -2^, 

K{ Ai 



where 

and uj = u>(Ki,K2) = —Kf — is the dispersion relation. 

Proof. It is sufficient to substitute the plane wave into the linear part of the matrix KP 
equation. ■ 



We use the transformation (|l,l| ) and introduce the following formal asymptotic Fourier 
expansion 

+00 

U(x,y,t) = £^ip n (£,ri,T;e)exp{i(nz)}, (2.1) 

n=— 00 

where z = K\x+K2y—ujt, j n = \n\ for n 7^ 0, and 70 = r is a non negative rational number 
which will be fixed later. The unknown functions ip n depend on e and it is supposed that 
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their limit for e — > exists and is finite; in the following this limit will be denoted with 
ijj n (^,r],T). Moreover we suppose that they can be expanded in power series of , i.e. 



MZ,V,r;e) =5>V£ ) (£,»7,t), Mt, V, r) = ^ (£, V, r). 
i=0 

We now introduce an analogous Fourier expansion 

+00 

W(x,y,t) = e 7n Vn(^r/,r;e)exp{i(nz)} 

n=— 00 

and obtain 



(2.2) 



( J FC 2 )(K 1 )-Vn + 0(e p ) 



In the following for simplicity we use the abbreviations tfi x = ^, ip { "( = <£, ip^' = A (and 
^ 0) = <f>n, = «)■ 

The final goal is to obtain the evolution equation satisfied by the modulation amplitudes 
Vl/ = ^f(^,r],T) and ^(^,tj,t) and to understand how it is modified by choosing different 
wave numbers. 



i,(0) 



,(0) 



Proposition 2.1. TTie matrix system (l.i )-(1.4) can be obtained applying the reduction 
technique to the matrix Kadomtsev-F 'etviashvili system (l.i). 



Proof. We insert the expansions (|2.l| ) and (|2.2| ) into the equation ( |l.2|) and consider the 
different equations obtained by considering the coefficients of the Fourier modes. 

It is convenient to separate the contributions of the linear and nonlinear parts by writing 

F 7n D ib — r 2 F 

where D n is a linear differential operator acting on ^ n (^,r/,T) and F n is the contribution 
of the nonlinear part. The operator D n is 

D n = {-inu + e q d T - V^d^ - V 2 e p d ri ) + (inK x + e p d^f - (inK 2 + e p d v ) 
- {i/K x ) (e p d v - {.K 2 /K x )d{) + {1/Kf) e 2p (d^ - {K 2 /K x )d^) . 

F n can be derived, by assessing the importance of the different terms, which originate 
from the nonlinear interaction of the Fourier amplitudes VVi(£> ^ r ) : 

F 2 = 6iK x ^ 2 + 0(e p ), 
F = e p(3{^,$} 5 -n/3[(/4 p) ,$] -iv^^S,*]) +0(e 2p ,e 2 ) , 

Fa = £ r ~ l (-iVS[n, ¥] - *V3^[*, A] + 3i^i{A, ^ 

+ 3ieK 1 {V'2,$} + 0(e r+p " 1 ,e 3 ) , 



and so on. 
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By setting q = 2,p = l,r = 2 for the proper balance of terms, we obtain the equations 
for the Fourier modes at the lowest order for n = 0, n = 1 and n = 2: 



ZKl - |§) A 5 + 2^A, - fi„ - 3{*, 



+ i\/3 



+ iv / 3 



* T + i 3iTi + 



#1 



= 0, 
2^* 



if. 



2 5^ 
1 



+ «V3[0, *] + iv^-^t*, A] - 3iKi{A, - 3iKi{^ 2 , $} = 0. 



Finally, after the cosmetic rescaling 

VSKxA -> A, \/3n 





d> A=^ 



??' = Kir/, 



(2.3) 



we arrive at the matrix system of nonlinear evolution equations (|l,3j)-( fL~4] ). 



This matrix system must be integrable by the spectral transform, because it has been 
derived from an S-integrable equation. This is explicitly demonstrated in the next section. 

Let us now look in more detail at the integrable NPDEs implied these results. If we 
take $ = ty* , N = 1, we obtain the equation 



i* T + Li* + X* = 0, L 2 x = 2Li|^r, 



(2.4) 



with 



Li = 
L 2 



Xv = 2|*g + 2V3% 
(3 + A 2 ) — + 2A 



d 2 



a 2 ^ 

d^drj dr] 2 ' 

d 2 d 2 



d^dr/ dr] 2 



The NPDE ( |2.4D , up to trivial rescalings, coincides with the Davey-Stewartson equa- 
tion [6], whose integrability is well known [7]. Note that the S-integrable equations found 
in [10-11] are different from the standard Davey-Stewartson equation and then not con- 
nected with the S-integrable system ( |l.3|) -(fL~4|). 

In the case N = 2, we get a nonlinear system for eight interacting fields. However, an 
interesting reduction is possible, if we set 



4>2 tpl 



<Pl <P2 
<P2 <Pl 



A 



Ai A 2 
A 2 Ai 



n 



n 2 sii 
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From the matrix system ( |1.3| )— (1.4), we obtain \P = $* and 

2 



iih,r + Inih + 2V3(A^! + A 2 ^ 2 ) - 2 (l^ilVi + + 2Vi|^2 
^ 2) T + L x ih + 2V3(Ai-0 2 + Aa^i) - 2 (|V> 2 | 2 ^2 + + 



0, 
0, 



L 2 Ai = L 3 (|Vi| 2 + |V> 2 | 2 ) , ^ 2 A 2 = L 3 (ViV4 + ^1) , (2-5) 



where 



L x = - (3 + A 2 ) <9f + 2A5|„ - a 2 , 
L 2 = (3 - A 2 ) df + 2Aa|„ - a 2 , L 3 = 2 v / 33f . 

Integrable Davey-Stewartson type equations and system of equations have been exten- 
sively investigated by many authors [16-20]. A detailed list of Davey-Stewartson systems 
and equations integrable by the inverse scattering method has been recently given [21]. 
The system of equations ( fEE| ) does not appear in these papers. We expect that this new 
system be integrable by the inverse scattering method, because it has been obtained from 
an integrable equation and the property of integrability is expected to be maintained 
through the application of the reduction method. The integrability of the system of equa- 
tions (|L~S|)- (|1.4| ), and of the system (|2.5D which is a particular case, is demonstrated in 
the next section. 

3 The Lax pair for the integrable system of equations 

In this section we apply the reduction method also to the Lax pair of the starting matrix 
equation ( |1.2j ) , to demonstrate explicitly the integrability by the spectral transform of the 



matrix system (1.3)— (1.4), and we thereby identify the Lax pair for the system of equations 



Lemma 3.1. The Lax operators (L,A) of the matrix KP equation are 

i cP' 

L = ^=— + ^-U(x,y,t), L<j)(x,y,t) = 0, (3.1) 
d 3 d 

A = 4 7n? - 6C/ ( X > y> *)«- - 3C/ *( X > y> *) + *>/3w; (3.2) 

ox s ox 

with 

4> t (x,y,t) + A<j)(x,y,t)=Q. (3.3) 
Proof. It can be verified by direct substitution that the operator relation 

In =i[L,A] = i(LA — AL) 
reproduces equation (|L~2j). ■ 



Proposition 3.1. The matrix system (l.S)-(1.4) is S-integrable and its Lax pair (L,A) 



is 



Lcf) = 0, (3.4) 
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where 



L= {l 2 1 1 L22 ) ' Lu = I ( id " + " x)d t) ' Ll2 = 

L21 = ^22 = I (ify - i(V3 + A)%) , 4> = (^ + _ ) , (3-5) 



/ is the N <S> N unit matrix and 



where 



A 12 = -2V3^a c - (Vs + a) * € + 

A^3 



vi 22 = -6»ia£ + in + 1 (y$ - a) a + % 



12 



1 



Proof. Let us apply the reduction method to the Lax pair (|3.1|) - (|3.3[) of equation (JL 

The components (j)j(x,y,t), j = 1,...,N, of the column vector <p(x,y,t) can be ex- 
panded in Fourier modes as follows 

+00 

<f>j(x,y,t) = ^2 e7 " ( Aj>(^ r /' r ; e ) ex P i + + ^t) + , (3.7) 

n=— 00 ~~ 

where z = -ftTix + K 2 y — Lot, the 0j iri (£, 77, r; e) depend parametrically on e and remain 
finite when e — ► 0, the 7„ are non negative rational numbers and A m , m = 1, ... ,3, are 
real constants to be determined. 

Inserting now the expression for 4>j(x, y,t) in ( |3.l[) , we derive a series of relations which 
are generated by the coefficients of the Fourier modes. Each relation must be valid for a 
given order of approximation in e. 

In particular, for the fundamental harmonics n = ±1, considering terms O in ( |3.1[) 
and (13.31), we obtain 



i ( iK 2 , A / iifi , A 



T^r + *A 3 +4 ±— 1 + zAi =0, 



and then 



A- K2 A- ^(XZ+kA A- K * ^KK 
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We thereby understand that the harmonics 

<t>j,i, j = l,...,N, (3.8) 

are fundamental, i.e. for them 7„ assumes the smallest value, 7„ = 0. 

The successive order e for the equation (|3.1|) allow us to obtain the new spectral prob- 



lem, because all the (pj tU may be expressed by means of the fundamental harmonics (3J 
which are connected through the relations: 

+ i(V3-\\ <f> +4 - = 0, (3.9) 
i(f>- r, - i (V3 + \) £-$<!>+ = 0, (3.10) 

where we set (c^i; j = 1, . . . , N) = <p+, (<pj,-i; j = 1, . . . , N) = </>-. 

By means of the variable rescaling Q2,3| ), and by introducing the 2N <8> 2N matrix 
operator L, we arrive at the final form (|3.4|)- (|3.5D . 

To calculate the temporal evolution, we must insert the expression Q3.7| ) in ( |3.3|) and 
consider the relation obtained for the different harmonics n and for a given order of 
approximation in e. If we consider the first order in e, we obtain again the spectral 
problem (Op— fl3T5|) . Only if we take into account the next orders of approximation of 
equation (3^), i.e. the order e 2 , the temporal evolution can be determined. However, new 
quantities, the corrections 0±(£, rj, r) of order e to the fundamental harmonics <j>±(£, r), r), 
appear. These unknown quantities can be eliminated in the equation (|3.3|) by taking 
advantage of the relation obtained from equation ( |3.1|), considering terms of order e 2 . 
This elimination is possible only because equations ( |3.l|) and (|3.3| ) are identical at the 
order e. In particular, if we consider ( |3.3| ) calculated to the order e 2 for n = ±1, we get 



+iT + 12i + A^ - 6*0_, f - 3* c 0_ + ^ (V, - 0_ 
+ iVsn - 6iA + Ai^ 4>+ + i f-6Ai - QK X + v^^J $<f>+3 

- — + 1 r 1 - wi) *+* ~ *+-) = °' 

<l>-,r + 12* (-^ + A^ - 6*0+^ - 3$^+ -^U v - 0+ 
+ iV3fi - 6iA ("^ + A i) 0- + * (" 6A i + 6i ^i + ) ^^- 3 

To evaluate this expression we took advantage of the fact that <j>±$ are connected with 
the fundamental harmonics (these relations are obtained from ( |3.lD for n = ±3 at the 
lower order in e): 



5+3 



gj) «-3 = (g, ) (3.11) 



The Matrix Kadomtsev-Petviashvili Equation 



19 



We now consider the equation (|3.l| ) at the order e 2 for n = ±1, which provides the 
corrections (j>+(£, rj, r), f], r). Via the transformation ( |2.3D and after a lengthy calcu- 

lation we arrive at the final form ( |3.6[ ) for the 2iV <S> 2N matrix operator A, which satisfies 
the equation 

4> T + A<i> = 0. (3.12) 



The determination of the Lax pair (|3.5|) and (|3.6[), which satisfies the equations (I 



and (3.12), demonstrates the S-integrability of the system (|1.3|)-(|L4|). 



4 Conclusion 

We have derived a new, integrable, and presumably of applicative interest, nonlinear mat- 
rix system of evolution equations of Davey-Stewartson type from the integrable matrix 
equation ( |1.2j) , by means of an extension of the reduction method based on Fourier ex- 
pansion and space-time rescalings. It reduces to the standard Davey-Stewartson equation 
in the single mode case and to a new integrable system of two interacting fields in the 
N = 2 case. Moreover, we have applied the reduction method to the Lax pair (|3.1|)- 
(|3.3| ) of the original equation and have demonstrated the integrability property of the 
new matrix system of equations, by exhibiting the corresponding Lax pair (|3.8| )-( |3!9|) and 

We have outlined the approach that permits to obtain such system of equations and 
the next steps will be the explicit resolution of the spectral problem and the possible 
identification of localized or asymptotically finite solutions. 

It is also convenient to push the approach beyond its "leading order" application by 
considering different rescalings in the two spatial variables or looking at special cases when 
some key parameters vanish, in analogy to the case of the model equations treated in [8-9]. 
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